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An approach is suggested defining effective sums of divergent series in tfie form of self-similar 

■ exponential approximants. The procedure of constructing these approximants from divergent series 
with arbitrary noninteger powers is developed. The basis of this construction is the self-similar 
approximation theory. Control functions governing the convergence of exponentially renormalized 
series are defined from stability and fixed-point conditions and from additional asymptotic conditions 
when the latter are available. The stability of the calculational procedure is checked by analyzing 

■ cascade multipliers. A number of physical examples for different statistical systems illustrate the 
O ' generality and high accuracy of the approach. 
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I. INTRODUCTION 



> 

B 

. Summation of divergent series is the problem of great importance in theoretical physics, apphed mathematics and 
O ■ engineering. This is because realistic problems are usually solved by means of some calculational algorithm often 
^ O^ . resulting in divergent sequences of approximations. Assigning a finite value to the limit of a divergent sequence 
is called renormalization or summation technique. The most widely used such technique is Pade summation [1]. 
However, the latter has several shortcomings. First of all, to reach a reasonable accuracy of Pade approximants, one 
^ ' needs to know at least in the magnitude of ten to twenty terms of a perturbation series. Unfortunately, so many terms 
are often not available because of the complexity of a considered problem. Second, Pade approximants are defined 
for the series of integer powers. But in many cases asymptotic series arise having noninteger powers. Third, there are 
quite simple examples [2] that are not Pade summable even for a sufficiently small variable. Moreover, the standard 
: Pade approximants do not converge at infinity, since infinity is an essential singularity [1] . The latter deficiency can be 
QQ ' sometimes treated with the help of two-point Pade approximants [3] . But these can only be constructed when there 
^\ are two perturbation expansions in the vicinity of two points, so that such expansions have compatible variables, which 
often is not the case, especially when one of or both these expansions contain noninteger powers [4]. For instance, only 
I rational powers can be described at infinity [3]. This is because a Pade approximant is a ratio of two polynomials, say, 
^ of m and n order. Therefore, when the variable tends to infinity, the asymptotic behaviour of this Pade approximant 
is of power law with the power m/n, which is a rational number. Fifth, in many physical problems the quantities 
of interest exhibit at infinity exponential behaviour, which in principle cannot be described by Pade approximants. 
One more well known difficulty when dealing with Pade approximants is the appearance of spurious poles [1]. The 
Q last, but not the least, Pade summation is rather a numerical technique providing answers in the form of numbers. 
Therefore, it is difficult, if possible, to analyze the results when the considered problem contains several parameters 
to be varied, since for each given set of parameters one has to repeat the whole procedure of constructing a table 
of Pade approximants and selecting of them one corresponding to a visible saturation of numerical values. Being 
^ I a numerical technique, Pade summation shares the difficulties of other numerical methods, like nonlinear sequence 
transformations, and sometimes is less effective than the latter [5,6]. 

The aim of the present paper is to develop an analytical method, free of the Pade approximation difficulties, for 
summing divergent series containing any number of terms ( just a few or many) with arbitrary noninteger powers. The 
method is based on the ideas of the self-similar approximation theory [7-14] in its algebraically invariant formulation 
[12-14]. The novelty of this paper is in what follows: 

(i) The approach is generalized by constructing self-similar exponential approximants from the series with arbitrary 
powers, integer as well as noninteger. (ii) It is shown how such exponential approximants can be made compatible 
with additional asymptotic and boundary conditions, (iii) Stability and fixed-point conditions are discussed and 
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concrete prescriptions, for defining control functions and for checking the stability of the calculational procedure, are 
formulated, (iv) It is demonstrated that the approach is applicable in all cases, when either just a few terms of a 
series are known or when many terms are available, (v) We emphasize that in all the cases, even when a large number 
of terms of a series are involved, we obtain analytical formulae that are convenient for considering with respect to 
a change of physical parameters. The possibility of deriving analytical expressions is a characteristic feature of our 
approach differing it from numerical methods. 

Also, we illustrated the generality of the approach by applying it to several interesting physical problems of quite 
different nature. Some of these problems could be treated by Pade summation or by the standard renormalization- 
group technique. When such treatments have been done, we compare their results with ours. In those cases for 
which numerical data are available, we estimate the accuracy and show that this accuracy is either comparable or 
higher than that of other more complicated resummation techniques. We construct as well self-similar exponential 
approximants for several physical problems for which other summation methods have not been applied because of 
technical difficulties. 



II. NONINTEGER POWERS 



The main idea [7] of the self-similar approximation theory is to put into correspondence to a sequence of perturbative 
terms a dynamical system for which the approximation number would play the role of discrete time. Then the 
transformation from one approximation to another can be represented as the evolution of this dynamical system [7] . 
The corresponding evolution equation may be formulated as the property of functional self-similarity, which is a 
necessary condition of fast convergence [8]. A dynamical system with discrete time is called a cascade. A fixed point 
of a cascade, representing a sequence of approximations, corresponds to an effective limit of this sequence [7,8]. To 
guarantee convergence, the fixed point must be stable and the sufficient conditions of stability can be formulated in 
terms of multipliers [9-11]. Additional possibilities open if we require that the procedure of finding an effective limit of 
a sequence is invariant with respect to algebraic transformations [12-14], which permits us to deal not with the initial 
sequence of approximations but with a sequence of its transforms. General ideas and the mathematical foundation 
of the self-similar approximation theory have been described in detail in our previous papers [7-14]. Not repeating 
them here, we begin at once with considering the case when for a sought function f{x) one derives an approximate 
perturbative expression 

k 

Pk{x) = '^an x"", (1) 

n=0 

in which are arbitrary real numbers, integer or noninteger, positive or negative, but with the sole requirement 
that they form an ordered sequence {a„}, that is, either strictly increasing or strictly decreasing sequence of terms. 
Following the method of the algebraic self-similar renormalization [12-14], we define the algebraic transform 

k 

Pk{x, s) = pk{x) = Y,an , (2) 

n=0 

where s is real. Then, by means of the equation 

Po{x,s)=aox'+''° (3) 



we obtain the expansion function 



X l/(s+ao) 

x(^,s)=(J) . (4) 



Substituting the latter into Eq. (2), we have 

yk{v,s) = Pk{x{ip,s),s) = '^an ( ) . (5) 



k f if^ is+a„)/{s+ao) 



n=0 



ao 



The family {yu} of transforms (5) is called [10, 11] the approximation cascade, since its trajectory {yi^{ip,s) \ k = 
0, 1, 2...} is bijective to the sequence {Pk{x, s) | fc = 0, 1, 2...} of approximations (2). A cascade is a dynamical system 
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in discrete time fc = 0, 1, 2..., whose trajectory points satisfy the semigroup property yk+p{(p, s) = yk{yp{(p, s), s). The 
latter equation is a particular type of functional equations [15] and in dynamical theory it is related to autonomous 
dynamical systems [16]. The physical meaning of the above semigroup relation can be understood as the property of 
functional self-similarity [17] with respect to the varying approximation number [7-11]. The self-similarity relation is 
a necessary condition for the fastest convergence criterion [8,9]. 

For the approximation cascade {yh}, defined by transform (5), the cascade velocity can be written as a finite 
difference 



Vkiv, s) = yk{<P,s) - Vk-iif, s) = ak 



(s+afc)/(s+ao) 



This is to be substituted into the evolution integral 

rPu 



/ 



dip 



r, 



(6) 



(7) 



in which Pk = Pk{x, s) and r is the minimal time needed for reaching a fixed point = -P^ (a;, s, r). Integral (7) with 
velocity (6) yields 



Pk{x,s,T) = 



Pk-ii^.s) 



V akT 



-l/v 



where 



u = Vk{s) = 



Taking the algebraic transform, inverse to Eq.(2), we find 

pI{x,s,t) = x~^Pkix,s,T) = 



ctk - (Xq 
s + ao 



-v ( \ _ ^ "'kT 
Pk-1\^) 1+1/ ^ 



(8) 



(9) 



(10) 



Exponential renormalization [13,14] corresponds to the limit s — > oo, at which 

lim ffe(s) = 0, lim svkis) = ak- ao- 

s—*oo s— *oo 

Then Eq. (10) gives 



lim pI{x,s,t) =pk-i{x)exp ( — rx"*" "° 



(11) 



Accomplishing exponential renormalization of all sums appearing in expressions of type (11), we follow the bootstrap 
procedure [14] according to the scheme 

Pk{x) ^pI{x,S,t) Fk{x,Tl,T2,...,Tk), (12) 

with k > 1. 

Let us illustrate explicitly how this exponential bootstrap works. The initial approximation from Eq.(l) is 

poix) = aox''° . (13) 

If we limit ourselves by the first-order term 

pi{x) = po{x) + aia;"S 
then, following the renormalization scheme (12), we get 

Fi(a;,Ti)=po(a;)exp(6ia;''^), (14) 

where 
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bi = —Ti, (3i = ai-ao. (15) 
ao 



Involving the second order term 

P2{x) = pi{x) + a2x"^, 

we find 

F2{x,Ti,T2) =po{x)ex.p {bixf^^ exp {b2X^^)) , (16) 

with the notation 

62 = — T2, 02 = 0:2- ai. (17) 

Continuing the same procedure, for the k—th order expression (1) wc obtain 

Fk{x,Ti,T2, ...,Tk) = po{x) exp [bix'^^ exp [b2x'^^ ... exp [bkX^'')) ...) , (18) 

where 

bk= — —Tk, I3k = ak-0ik-i- (19) 

The quantities r„, with n = 1, 2, . . . fc, in the rcnormalizcd form (18) play the role of control functions [7-14]. In 
the following section we shall show the ways of defining these control functions. Assume, for a while, that the latter 
have already been defined giving r„ = t„(x). Substituting these functions into Eq.(18), we come to the self-similar 
exponential approximant 

/fc*(x) =^^fc(x,Tl(x),T2(x),...,Tfe(,T)). (20) 

Note that the constructed approximants (20) are different from the iterated exponentials introduced by Euler [18]; 
the properties of such iterated exponentials being reviewed, e.g., in Refs. [19,20]. 

Although in this paper wc shall deal with physical problems related to scries with powers being real numbers, 
nothing prevents us from generalizing the whole approach to series with complex powers. Such complex exponents 
appear in the problems with discrete scale invariance [21,22], which has recently been documented in the models of 
rupture, earthquake processes, financial crashes, in the fractal geometry of growth processes, and in several random 
systems [21-24]. Our approach can be straightforwardly applied to series with complex powers. Then, the variable s 
in the algebraic transform (2) has also to be considered as complex. Therefore, the corresponding control functions 
Sk{x) become complex. However, before passing to the general case of complex powers and, respectively, of complex 
control functions, we need first to develop the approach for arbitrary real powers and real control functions. 



III. CONTROL FUNCTIONS 



The role of control functions, by their definition [26], is to provide convergence for initially divergent sequences. 
There are several ways of incorporating control functions into an iterative algorithm and of defining them. One way 
of introducing these functions is by including them into an initial approximation. The most often used definitions of 
control functions are through the minimal-difference [26-28] or minimal-sensitivity [29-36] conditions. In some very 
simple cases, like zero-dimensional and one-dimensional oscillators, when high-order terms of perturbation theory are 
available explicitly, control functions can be found from the direct observation of convergence of this theory [37-40]. 
All these variants are particular types of quasifixed-point conditions [7-11]. 

Another way of introducing control functions is through an algebraic transformation [12-14]. In addition, the 
minimal time appearing in the evohition integral (7) under each rcnormalization step plays also the role of a control 
function. After k steps of the rcnormalization procedure, the self similar approximation (20) contains k time like 
control functions r„, n = 1, 2, . . . fc. 

The simplest way of defining these control functions would be to remember that the effective time in integral (7) 
corresponds to the minimal number of steps needed for reaching a fixed point. When no other restrictions are imposed, 
the minimal number of steps is, clearly, equal to one. Putting r„ = 1 for all n = 1, 2, . . . fc in Eq.(20) gives 

rk{x)=Fk{x,l,l,...,l). (21) 
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A more elaborate definition of the time-like control functions can be formulated by involving one of the variants of 
fixed-point conditions. To this end, let us put r„ = 1 for all n = 1, 2, . . . fc — 1, except the last step for which t„ = r 
is yet undefined. Consider the sequence {fk} consisting of the terms 

fk{x,T)=Fk{x,l,l,...,T). (22) 

Following the standard procedure [10-14], it is possible to construct an approximation cascade with a trajectory 
bijective to the sequence {fk}- For the cascade velocity we may write the finite difference 

Vk{x,T) = fk{x,T) - fk-l{x,T). (23) 

When approaching a fixed point, the cascade velocity tends to zero. Therefore, the condition to be as close to the 
fixed point as possible is the minimum of velocity (23). This condition 

min\Vk{x,T)\ = \Vk{x,Tk{x))\ (24) 

T 

provides us the definition of the time-control function Tk{x). Eq.(24) is the general form of the minimal- velocity 
condition [12,41]. In particular, taking into account definition (23), we may have the equation 

fk{x,T)-fk-i{x,T)=0, (25) 

whose solution r = Tk{x) is the sought control function. For the exponential approximants (18), Eq. (25) yields the 
equation 

r = eicp ( xl'A . (26) 
\ak-i J 

The solution to Eq. (26), that is, r = Tk{x), being substituted into (22), leads us to the self-similar approximation 

mx) = fk{x,Tkix)). (27) 

The described scheme of defining control functions is applicable when no additional restrictions are imposed on the 
behavior of the sought function /(x). It may happen, that, in addition to expansion (1), an asymptotic behavior of 
f{x), as X ^ xo, is known. Then the asymptotic condition 

fix) ~ fo{x), x^xo (28) 

plays the role of an imposed restriction. And the constructed self-similar approximations are assumed to satisfy the 
asymptotic condition (28). In such circumstances, some of control functions are to be chosen so that condition (28) 
be valid. This can be done in the following way. 

Let us renormalizc a scries pk-i{x) to a self similar approximation f^_i{x) with time-like control functions defined 
according to a scheme described above. Limiting ourselves by such a (A: — l)-step renormalization, we obtain from 
Eq.(lO) 



Fk{x,s,T) 



"0 



(29) 



with the same notation (9). Now we require that the obtained expression (29) would satisfy the asymptotic condition 

-Ffe(a;,s,T) ~ /o(a;), x ^ xq, (30) 

in accordance with (28). The control functions s = Sk{x) and t = Tk{x) are to be chosen so that condition (30) be 
valid. With the so defined control functions, we come to the self-similar approximation 

f^ix)=F^{x,Sk{x),Tk{x)) (31) 

possessing the desired asymptotic property (28). 

To clearly illustrate the latter variant of defining control functions, consider the integral 

1 f°° 

J{g) = -j= J exp (-x^ - gx^) dx, (32) 
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which has the meaning of the partition function for the zero-dimensional anharmonic model, where g is called a 
coupling parameter. The series for integral (32) are frequently used as a model for strongly divergent perturbation 
expansions in quantum field theory [25]. 
The weak-coupling expansion of (32) gives 

J{g)c^a + bg+..., {g ^ 0), (33) 

with a = l, 6=— |. In the strong-coupling limit one has 

J{g)^Ag-^/'' + Bg-^/'' + Cg-'/' + Dg-y\.., (<? - oo), (34) 

where 

1.813 0.612 0.227 
A= — B = 6 = — 

\'TT ^/tT a/tT 

The strong-coupling expansion (34) is of the form of series (1) with the coefficients 

ao = A, ai = B, 02 = C, 03 = -D 

and with noninteger powers 

13 5 7 

ao^~-, ^1 = —^, "2 = -4' "3 = -^. 

We shall renormalize expansion (34) subject to the limiting condition 

lim J{g) = a, (35) 

which follows from (33). 
Starting from 

Mg) = Ag-^l\ (36) 

we get 

^2 (5) = U9) exp (l^"'/' exp (§5"'/')) • (37) 
At the next step, according to (29), we have 

F3*(ff, s, r) = Vr^{g)r + 7 r\ '"'^ , (38) 



with 



6 V Dt 



l-4s' ' 
In the weak-coupling limit, Eq. (37) yields 

j;{g) Ag-'/^ (3^0) 

and Eq. (38) gives 

F^{g,s,T) = {A--g"'^ + ^g^''y"\ 
The latter expression satisfies the limiting condition 



lim F*{g,s,T)= a, (39) 
9^0 
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corresponding to (35), if and only if s = and 7 = a or, respectively. 



v = 6, 



-1/6 



Substituting the found s and r into (38), we get 

^3(5)= [(^2(.9))"' + a-6' 
Using here form (37) and remembering that a = 1, we finally obtain the self-similar approximation 



^3(5) 



1 



^3/2 



/~6B 



exp 



C 



-1/6 



(40) 



corresponding to (31). 

The exponential approximant (40) represents integral (32) with a very good accuracy. The maximal percentage 
error is 2.7% occurring at g = 0.1. And the error for the physically most interesting region of intermediate g Ki 1 is 
only 0.7%. 

Consider another simple example evidently illustrating the generality of our method that works when other methods 
do not. Assume that we need to find an approximation for a bounded function f{x), with — 00 < a; < 00, having the 
following asymptotic properties: 

/(a;) ~ 1 + a; + 0(a;2) (a; ^ 0), 
f{x) ~ e"^ (a; ^ -00), 



■^(^)-^' (a; ^+00), 

where a is an irrational power. To our understanding, Padc technique is principally inappropriate in this case. While 
in our method, accomplishing the same procedure as is explained above, we easily obtain 









exp ^ 













-a/2 



IV. STABILITY CONDITIONS 



In order to check the stability of calculational procedure, one has to analyze mapping multipliers [9-11]. Several 
kinds of multipliers occur in the process of construction of self-similar approximations, each kind being related to the 
corresponding approximation cascade. 

The first approximation cascade, appearing in our investigation is the cascade {uk} composed of transforms (5). 
The local multipliers for this cascade are defined as 



9V {s + ao)aQ \aoJ 

The image of multiplier (41) in the a;-space can be obtained with the use of relation (3), which yields 



(41) 



mk{x,s) = iJ,k{Po{x,s),s) 



dPk{-< 



\/()x 



dPo{x, s)/dx 



^ (s + ao) ao 



(42) 



The trajectory {yk{v, s)} of an approximation cascade {yk} is locally stable [11] at the fc-step if \nk{f, s)\ < 1. Note 
that the case of |/ife| = 1 is called neutrally stable and that of = can be termed superstable. If \^k\ < 1 for 
some s and all ip from a given domain, then |mfc(a;, s)| < 1 for the same s and all x from a domain defined by relation 
(4). When the trajectory of a cascade {yk} is locally stable for all fc = 0, 1, 2, . . ., that is, when \mk{x, s)| < 1 for all 
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fc, then the sequence {P/j(.t, s)} converges uniformly with respect to x. The local stability for all fc = 0, 1, 2, . . . can 
be called the global stability. The global stability of a trajectory is a sufficient condition for the convergence of the 
corresponding sequence [11]. But it is not a necessary condition. Thus, the sequence {Pfe(a;, s)} may be convergent, 
but the cascade trajectory not everywhere locally stable, so that the stability condition \mk{x, s)\ < 1 becomes valid 
for all k starting from some ko, but for k < ko this condition may be broken for some k. 

The concept of stability suggests a recipe of defining the control function s = Sk{x). The latter can be defined so 
that to minimize the absolute value of multiplier (42), which can be named the principle of maximal stability [12-14]. 
Substituting the found Sk{x) into Pk{x,s), we obtain Pk{x,Sk{x)). The renormalized sequence {Pk{x, Sk{x))} can 
become convergent even if the initial sequence {Pfe(.T, s)} was not. Defining local multipliers for the new sequence, it 
is necessary to take into account that Pk{x, Sk{x)) depends on x explicitly as well as through Sk{x). The corresponding 
multipliers are 

, . _ dPk/dx + jdPk/dsk) {dsk/dx) 

dPo/dx+{dPo/dsk){dsk/dx) ' ^ ""^ 

where Pi~ = Pk{x, Sfc), Sfe = Sk{x), and the partial derivatives mean that another variable is kept fixed. In general, 
mk{x) differs from (42). But if Sfe(x) is a slowly varying function of x, such that the derivative dsk/dx in mk{x) 
can be neglected, then mk{x) approximately coincides with mk(x, Sk(x)). The stability condition |mfe(a;)| < 1 for all 

= 0, 1, 2, . . . implies the convergence of the sequence {Pk{x, Sk{x))}. Let us stress again that stability is a sufficient 
condition for convergence, but not necessary. The sequence {Pk{x, Sfe(.x))} may converge even if the stability condition 
is not valid for a finite value of k. Moreover, the convergence of the sequence {Pk{x, Sk{x))} is not compulsory for 
us. This sequence is not yet the final product of the procedure but is to undergone the dynamical renormalization 
involving the evolution integral (7). 

After the multiple dynamical renormalization, in the exponential variant we consider here, we come to the sequence 
{Ffe(.T,T)} of terms given by expression (18), with the shorthand notation r = {n, r2,..., Tfe}. The stability of the 
corresponding trajectory is characterized by the multipliers 

M.(.,r). ^!:f'-;/f . (44) 
dFi{x,T)/dx 

The stability condition \Mk{x,T)\ < 1 for all A: = 1, 2, . . . guarantees the convergence of the sequence {Fk{x,T)}. 
This tells us that the time-control functions Tk{x) could be chosen so that to minimize \Mk{x,T)\ . In particular, the 
minimal value of (44) can be zero. Then the equality Mk{x,T) = 0, under the assumption that dF\{x,T)/dx ^ 0, 
yields the equation 

—Fk{x,r)=0 (45) 

restricting the choice of time-control functions. 

Defining all time-control functions by one of the ways discussed above, we obtain, as a final result, the exponential 
approximant (20). The sequence {fk{x)}, with k = 1,2, of these exponential approximants is what we ultimately 
need to analyze with respect to its convergence. For that purpose we can formulate sufficient conditions for convergence 
studying the stability of the corresponding cascade. To this end, following the standard procedure [7-11] , we define 
a function x{(p) by the equation 



/*(x) = <^, x = x{^). (46) 



Then we introduce the transformation 



yl{^) = f*k{x{^)). (47) 
The family {yl} of the transformations introduced in (47) forms a cascade. The local multipliers are given by 



M^(^) ^ ^yU^). (48) 



The image of (48) in the a;-representation is 



_ df;:ix)/dx 
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The cascade trajectory is locally stable at a fc-step when 

\ml{x)\ < 1. (50) 

If this condition is valid for all k, then the sequence {/^(a;)} converges. 

The validity of the stability condition (50) is what finally justifies the whole rcnormalization procedure guaranteeing 
the convergence of the renormalized sequence {f^{x)}. This holds true irrespectively of whether other multipliers, 
as (42), (43), or (44), satisfy the same condition (50) as the multiplier (49). The auxiliary multipliers (42)-(44) have 
appeared at intermediate stages of our multistep renormalization procedure. The minimization of these multipliers 
provides a recipe for defining control functions. This minimization is equivalent to the stabilization of local parts of 
a cascade trajectory. Remind that we start with series (1) which, in general, is strongly divergent. Accomplishing 
several stages of renormalization we step by step improve convergence properties, by making the related trajectory 
more and more stabilized. However, it is not necessary to demand that the trajectory would become completely stable 
at some intermediate step of the multistep procedure. The main is that the finally resulting sequence {/^(a;)} be 
convergent, for which it is sufficient to require the validity of condition (50). 



V. GROUND STATE PROPERTIES 



Now we pass to the consideration of physical examples illustrating explicitly how the method works. 



A. One-Dimensional Bose System 



The ground-state energy of the one-dimensional Bose system with the ^-functional repulsive interaction potential is 

known in a numerical form from the Lieb-Liniger exact solution [42] . We derive below a compact analytical expression 
for the ground-state energy e{g) as a function of the interaction strength g, valid for arbitrary g. In the weak-coupling 
and strong-coupling limits the following expansions are known (see e.g. Ref.[7] and references therein) : 

e{g) ~ + bg""/^ + eg"" + dg^^^^ + ... , {g ^ 0) , (51) 



a = 1, 



4 

3^' 



0.0654, 



-0.0018, 



while in the strong-coupling limit an exact result is available [43]: 

^2 



e{g) c^A= — 



(g^oo). 



(52) 



The expression satisiying both known limits, can be derived similarly to the example studied in Section III, except 
that in this case we start from the weak-coupling limit. The resulting approximant is 



el{g) = ag 



exp 



1/2 

-g ' exp 
a 



-3/2 



-3/2 



3/2 



-2/3 



(53) 



This expression works reasonably well in the region of g € [0, 10], where the exact numerical solution of the Bethe 
ansatz equations is available. The maximal error here is about 5%. 



B. Asymmetric Anderson Hamiltonian 

The single-orbital Anderson Hamiltonian describes the system consisting of localized d-electrons interacting with 
conducting s— electrons via a quantum mechanical exchange mechanism, whose strength is measured by V, which is 
the transfer integral between the s- and d-states. The energy U describing the Coulomb interactions between two 
rf-electrons is anottuT relevant physical parameter. The localized d— electron number rid is expressed through the 
self-energy E of d electrons at the Fermi level: 
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rid = 



1 1 



tanh " 



(54) 



where p stands for the density of states at the Fermi level for conducting electrons. In the case when d-level is fixed 
to the Fermi-level, the following expansion for S, in powers of the small parameter u = U/nA, |u| <C 1 , was obtained 
in Ref. [44]: 



E 
A 



-u (l + aiu + a2U^ + a3U^ + a^u^ + a^u^) (u — > 0), 



(55) 



ai 



-1, 



02 = 0.5326, 



03 = 0.6269, 



04 



-1.8071, 



05 = 1.027. 



A direct application of expansion (55), gives physically meaningful results in the region of —0.5 < u < 0.5, as is 
shown in Fig. 2 of Ref. [44]. However, the exactly known limits 



(u — > oo); 



rid 



1 (u ^ —oo), 



are violated when expansion (55) is used for the calculation of the impurity level occupancy (54). To our knowledge, 
the Padc approximants were not applied for resumming expansion (55), and it is not clear whether they can be applied 
at large values of the parameter u (see [2,45]). We use below the technique of exponential approximants and observe 
that they spontaneously recover the known limits. The following sequence of self-similar renormalized expressions 
can be written: 



; (exp(aiu)) , 



(56) 



exp I oiM exp I —u 

Oi 



(57) 



A =2^ 



02 / O3 

exp oiuexp — uexp — u 
' oi V 02 



(58) 



TT 

2" 



02 / 03 / 04 

exp oiuexp — uexp — uexp — u 
oi V V 03 



(59) 



TT 

— U 

2 



02 / 03 / 04 / 05 

exp oiuexp — ucxp — ucxp — licxp — u 
oi \a2 \a3 V04 



We observe, that already the third-order approximant {^)'^ leads to the result 

(ricda 



11 , -1 
tann 

2 TT 



(60) 



(61) 



which possesses the correct limits. The same is true when the higher order approximants 



{ricdl = 7;-- tanh ^ 



{rid)5 



11 1-1 
tanh 

2 TT 



are considered. An aposteriori analysis of multipliers suggests that the third-order approximant corresponds to the 

most stable trajectory. This could be expected beforehand, since the coefficients 04 and 05 are larger than one. 

The density of states of d-electrons as well as the resistivity R at zero temperature, could be also reconstructed for 
arbitrary u, e.g. 



(R{u)y 

\R{0)J 



H -1 



The shape of the curve {nd{u))'^ is very much like the smeared Fcrmi-distribution, while ' Einl 



-jKqj J has an asymmetric 

bell-shape. All artifacts, which appear at the curves, when the perturbative expansions are naively extended beyond 
the region of small u, are smeared out, and the renormalized curves appear to be rather smooth. 
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C. t— expansion 



The so-called expansion is a tool for a systematic improvement of variational calculations for Hamiltonian systems 
[46]. Using a t— dependent variational wave function and, after performing calculations, finally, taking the t ^ oo 
limit, one can hope to increase the quality of the variational estimate for the ground-state energy E. This idea is 
most frequently used in conjunction with heavy numerical calculations and Pade-approximants technique. In the case 
of the Id Heisenberg antiferromagnet, the expansion in powers of the parameter t was obtained explicitly [46], 

Ec^ -j-t + 2t^ + ^t^ -m'^ (t^O), (62) 

and the task of getting an estimate for the ground state energy from the asymptotic expression valid at t ^ 0, can 
be approached by the methods of Sections II-IV. Since the number of terms available is finite, it is not necessary to 
take the limit t ^ oo explicitly, but instead we can try to minimize the error caused by this inevitable truncation by 
demanding the minimal sensitivity of the renormalized expression E* with respect to the "time" t determining the 
"duration" of motion to the ground state energy, 

Applying the self-similar bootstrap in its super-exponential form, one can see that the solution to this equation does 
exist for any number of terms from the initial expansion (62). From the aposteriori analysis of the multipliers, we 
conclude that the most reliable value, corresponding to that obtained along the most stable trajectory, is 

E* = ~ exp ^At exp (^-2t exp ^ ^ ' (64) 

with t = 0.33 and E* = —0.446. Prom the viewpoint of an apriori analysis, it is admissible also to construct another 
seqiicncc of exponentials, not including the constant, mean field result. In this case, the answer for the ground state 
energy is —0.434. We conclude that moving along the two different, but stable, trajectories, we can determine E 
reliably, as E = —0.44 ± 0.006. The known Hulthen [47] exact result —0.4431 is located within these boundaries. 

VI. EFFECTIVE COUPLING 

A. Beta- Function of SU(2) Lattice Gauge Model 

The Callan Symanzik /?— function of the (3+1) dimensional SU(2) lattice gauge model in its weak coupling asymp- 
totically free regime may be presented in the form of an expansion in powers of the parameter g, where g stands for 
the coupling [48]: 

_M^ag2^6/ + ..., (g^O), (65) 
11 . 17 

a = 



247r2 ' 647r4 " 

In its strong-coupling limit, (3 can be presented as follows [48]: 

-^c^A + Bx'^ + Cx^ + Dx^ + Fx^, x=^ (5^00), (66) 

75' 625' 140625' 390625 

We will add to the expansion (66) one more trial term ~ Gx^'^ , and determine G from the boimdary condition following 
from the weak-coupling limit (65). Two starting terms from (66) will be renormalized to the form of an exponential 
approximant, mimicking an instanton contribution bridging these two limits [49], that is, a nonperturbative physical 
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mechanism coming into play in the crossover region, from strong-to- weak-coupling limit [50]. This term should 
disappear completely in the weak-coupling limit, guarantying rather sharp crossover. The last four terms will be 
renormalized to the form satisfying the boundary condition (65), in a way leading to a smooth matching of two 
limiting kinds of behavior. We obtain the following analytical form for the renormalized /3— function: 



, ^exp^-.^J+C.^^l-^-.^J .^J , (67) 

The shape of the resulting function is similar to those obtained in Refs. [48-50] . Crossover occurs at (/ ~ 1.2. Let us 
remind that the quality of the strong-coupling expansions is jeopardized by the interfering roughening transition [48]. 
We believe that the method discussed above, allows us to bypass this difficulty in a natural way. 



B. Kondo Effect 



We consider below an application of the exponential approximants to such an interesting problem as the Kondo 
effect [51], comparing our results with those of the field-theoretical renormalization group [52-54]. 

The behavior of the system, consisting of a local-impurity spin and conduction electrons, interacting by means of 
an antiferromagnetic exchange of strength J, changes from asymptotically free at high temperatures to that when the 
impurity is screened by electronic lump at low temperatures, via the crossover region whose onset is characterized by 
the Kondo temperature estimated as 



Tk = Dexp{- 



1 

2J 



(68) 



where D stands for the Fermi-energy of electrons. We consider below only the case of a single-channel Kondo model. 
Most of our knowledge about the problem came from the exact Bethe ansatz solution [55,56], from the field-theoretical 
renormalization group [52-54] , and from the Wilson numerical renormalization group [57] . 

Within the framework of the field-theoretical RG in its application to the Kondo crossover, the central role is 
played by the so-called invariant charge or effective electron-electron coupling J^v [52-54] , measuring the intensity of 
electron-electron interactions via the impurity spin: 



J 



l-F2Jln( ^ 

uj\ 



(69) 



here w stands for the typical external parameter of the problem, such as temperature, or magnetic field. If only the 
starting two terms from (69) are taken into account, the field-theoretical approach, through calculations based on the 
Gell-Mann-Low /3— function. 



(3 ~ -2 J2 



{J « 1), 



leads to the formally divergent, at T = Tk, expression for the invariant charge [52-54], 

J 



Jinv — 



1 - 2Jln 



(h) 



(70) 



(71) 



We apply below the technique of algebraic self-similar renormalization directly to the series (69) for Jinv, continuing 
them from the region of J <C 1 to the region of J ~ 1. When only two starting terms from (69) are taken into account, 
the optimal, from the viewpoint of stability conditions, solution is the following exponential approximant. 



7* 



Jexp 



2Jln( ^ 



■'In 



2J 



(72) 
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If we perform the self-similar renormalization with the control fimction s = 0, along the non- optimal trajectory, we 
will recover the expression (71). So, the fictitious pole is absent in our solution to the problem, although the typical 
energy scale determined from the condition 



2Jln( ^ 



1, 



coincides with the Kondo temperature (68). Expression (72) is formally divergent as w — > 0, in agreement with the 
conclusion of the numerical renormalization group [57,58] . 

More complicated situation arises in (quasi)-two-dimensional metal, when the Van Hove logarithmic singularity 
in the electron density of states can influence the Kondo effect. In this case, as was shown in Ref.[59], the full 
electron-impurity scattering amplitude F can be estimated as follows: 



J 



1 + Co J ( In ( 



+ ... 



Co > 0, (J < 1) , 



(73) 



i.e. the usual Kondo logarithm should be replaced by the squared logarithm, originating from the Van Hove singularity. 
The self-similar ly renormalized expression in this case again corresponds to the exponential approximant: 



J exp 



Co J In 



(74) 



and the characteristic energy scale, Kondo temperature, can be found from the condition 

2 



Co J I In 



leading to the estimate 



Tk = D exp 



const 



^/J J 



). 



1, 



const ' 



(75) 



Such a dependence of Tk on J was the main result of Ref. [59] , obtained as an outcome of a cmiibersome first order 
parquet summation . To our knowledge, the field-theoretical KG approach was not applied to the Kondo effect with 
the Van Hove singularity. On the other hand, the Bethe ansatz fails for this problem [59] . 

The higher order corrections were not considered in Ref. [59], because of the technical difficulties arising in more 
sophisticated parquet approximations. Our approach may be of interest in this context, allowing to find the corrections 
to Kondo temperature due to higher-order terms. Taking into account the higher-order perturbative terms, one has 



r = J 



1 + CoJ In 



+ cij2 In 



co>0, (J«;i), 



(76) 



with ci < 0. The following exponential approximant is optimal from the viewpoint of stability: 

2 " 



r = J exp 



Co J In 



exp(-Mj 



(77) 



Prom the condition 



we obtain the Kondo temperature: 



Co J In 



exp ( ^- J 

Co 



Tk = D exp 



cyvv2 



(78) 



This estimate suggests a decrease of Tk due to higher-order corrections. 
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VII. EQUATION OF STATE 



A. Classical Hard Spheres. 

The exponential approximants can be used for constructing equations of state for simple liquids. For the model 
system of hard spheres with the diameter d, widely used as a reference system, an empirical equation of state, suggested 
by Carnagan and Starling (see [60]), is 



nkT (1 - p)3 

connecting pressure p, temperature T, the number density n, and the reduced density p = irncP /6. It agrees very well 
with the molecular dynamics and virial expansion [60,61]. The theoretical virial formula, according to Percus-Yevick 
[60,61], is given as follows: 

p _l + p + p^-3p3_ ^^^^ 



nkT (1 - /?)3 

These two expressions almost coincide at low densities, e.g. at p = 0.1 the percentage error of Eq.(80), as compared to 
(79), equals —0.18%; while for the intermediate and high densities the agreement becomes very poor, e.g. at p = 0.5 
the percentage error is —15.385%, and at p = 0.7 it equals —37.141% . 
Consider the regular part of (80), defined as r : 

r~l + p + p2_3p3 (p^O), (81) 

as an asymptotic, low-density expansion for the regular part r(p), and try to continue the expression (81) from the 
region of p <C 1, to the region of p < 1. In order to extend the validity of (81), let us add to it one more trial term 
~ p^. It is reasonable to use for renormalization the last four terms from thus extended expansion for r, since the 
constant term describes the ideal gas behavior and we are interested in the region of high densities. Following the 
standard prescriptions of Sec.II-III, we write down the two exponential approximants, justified from the viewpoint of 
stability for the sequence of an aposteriori multipliers: 

4 (p, r) = 1 + p exp [p exp (-3pr)] , (82) 



rl (p, r) = 1 + p exp p exp (^-Sp exp ^p) ) . (83) 



We retained in expressions (82) and (83) the effective time r, introduced at the last step of the bootstrap procedure. It 
will work now as a control function r — t{p) determined from the minimal velocity condition mini- \r\ (p, t) — r'^ (p, t) | . 
For the sake of simplicity we choose a single control parameter r, instead of the control function, from the minimal 
difference condition imposed at a single point, chosen from the region of intermediate densities, say for p = 0.6. Then, 
r = 0.845. Recalculating 

P* _ rlip,r) ^g^^ 



nkT (1 - p)3 ' 

and comparing it with empirical formula (79), we obtain that at p = 0.1 the percentage error equals —0.112% ; at 
p = 0.5 the percentage error is —3.215%, and at p = 0.7 it equals —2.92%. Formula (84) is more accurate at high 
densities than our previous result [14], corresponding to 

P* _ rlipA) ^35^ 



nkT (1 - p)3 ' 

which works with the percentage error of —4.567% at p = 0.7. What is even more important here, is the possibility of a 
self consistent improvement of the quality of the equation of state, based on stability conditions. We should recall here, 
that phenomenological exponential-type expressions are well known in the theory of equations of state, beginning, 
probably, from the Hudleston equation (see [61] and references therein) and ending with its modern modifications 
[62-64]. For example, the Shinomoto's equation for the system of hard spheres reads: 



14 



nkT 



exp 



4p ( 1 + -p 



(86) 



and gives at p = 0.1 the percentage error 0.046%; at p = 0.5 the error is —6.289%, and at p = 0.7 it becomes 
-35.948%. 

A single-exponential approximation can also be obtained from Eq.(81) by our method: 



p* pexp (p - 3/9^) + 1 
^ " (1 - p)3 ' 



(87) 



giving at /9 = 0.1 the error -0.158%, at p = 0.5 the error is -14.498%, and at p = 0.7 it becomes -28.31%. 

We conclude that the multi-exponential formula (84) agrees well with the empirical formula (79), being superior to 
all other formulae in the region of high densities. 



B. Quantum Hard Spheres 



At low density p, the energy E for a boson system of N-h&vd spheres with the diameter c and mass m is known 
[65]: 



E ^2tt 
N ~ m 



c(l + Ci(pc3)'/' + ...), Ci 



128 

15V^ 



(p-0). 



(88) 



And as p ^ po, where po = V^/c^ is the maximal density for a system of hard spheres, the following expression is 
available [65]: 



N 2m 



A = 



2V3 



(P Po) : 



corresponding to a second order pole in the ground-state energy. Let us, in analogy with the previous example, 
extract the singularity, as p — > po, rewriting (88) as follows: 



E 2ttcP 



N m 



1/3 



1 + C^{pc'f' + ... 



(p-V3-po-^/^)' 



and, keeping only a few starting terms 



E 2nc 



„l/3 



N m 



1/3 



l-(^ 

Po 



1/3 



1/3- 



(p-0). 



(90) 



l-2(i-) +C,(pc^)^/V... 



(91) 



After the standard self-similar renormalization, involving two terms from (91), we obtain 



E* 27rc 



.1/3 



■ exp 



-2r 



Po 



1/3- 



(92) 



where the control parameter r should be determined from the known asymptotic formula (89), as p — > po- Finally, 



E* 27rc 



.1/3 



N m ( I/O -1/3 



■ exp 



1/3- 



1 1 M -1/3 
T = Pn 

2 KAttc^ 



(93) 



or, equivalently. 



E* 2-Kc 



N m 



,1/3 



(.-. 



/3-p-/3^H4.c- ; 



(94) 



Equations (93) and (94) should be compared with the empirical London equation of state [66]. Both equations give 
the results very close to each other. So, our derivation can serve as a justification for empirical formulae used for the 
system of Bose hard spheres. 



15 



C. Polymer Coil 



The expansion factor a of a polymer coil is represented as a function o? = o?{z) of the excluded volume variable 

z = B\fN (3/27rP)'^^^ , where N is the number of bonds of the length / each and B is an effective binary cluster 
integral [67-70]. The case of a polymer coil corresponds to ^ > 0. In the region of ^ -C 1, the perturbation theory in 
powers of z can be developed, giving the expansion [68] 



o?{z) ~ 1 + /ciz + k-iz^ + k-iZ^ + fc4z'* + fcsz^ + k^z^ {z 0), 



(95) 



-2.075385396, fcg = 6.296879676, 



ki = -25.05725072, fcj = 116.134785, ke = -594.71663. 
On the other hand, in the limit of z 3> 1, a is related to z by a simple power law 

a'^c^Kz'', 6=2(2i/-l) ^ oo), 



(96) 



where i> is the critical index, and K stands for the critical amplitude. One of the popular problems in the physics of 
polymer coils consists in the continuation of the expansion (95) to the region of arbitrary z. We derive below, using 
the self-similar renormalization, a simple equation of state for the polymer coil, valid for arbitrary z, and satisfying by 
design both known limits, (95) and (96). The coefficients in expansion (95), starting from fca, grow rapidly, so do the 
local multipliers. Because of this, we use only three starting terms from (95), stabilizing the renormalized expression 
by imposing the asymptotic condition (96). Finally, after the standard transformations analogous to those of Section 
III, we obtain 



{a\z)y 



exp ( kizexp [ —z 



3/b 



6/3 



(97) 



It is known, from different approaches [12,69,70], that 1/2 <z/ < 0.6 and 1.53 <K < 1.75. We take the values 
I' = 0.599 and K = 1.62, which we calculated in Ref.[12], and compare the equation of state (97) with the empirical 
equations of state [70] of Barrett-Domb, 



a 



. 20 , 2 

1 + —Z + AlTZ^ 
o 



1/5 



and of Yamakawa-Tanaka, 

= 0.572 + 0.428 (1 + 6.23^)^^^ . 
In the region < 2; < 12, our equation (97) yields a curve lying between those of eqs.(98) and (99). 



(98) 



(99) 



D. Ising Model 

The low-temperature expansion for the order parameter (magnetization) M of the three-dimensional Ising model 
on a f.c.c. lattice reads [71]: 



M{T)~l + hvy + cu'\ w = exp - 



(T ^ 0), 



(100) 



6 = -2, 



c = -24. 



Low temperature expansions arc non universal in the sense that they depend on the type of a lattice, spin etc. . .. At 
the critical point Tg, order parameter demonstrates a universal behavior, independent on the type of a lattice, spin 
etc. . .: 
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(101) 



where /3 « 0.325 is the critical index [72]. For the f.c.c. lattice, Tc ~ 9.8 [72]. We continue the low-temperature 
expansion (100) self-similarly to the whole region of <T < Tc, imposing relation (101) as a boundary condition. 
Following the standard prescriptions of Section III, we obtain 



M*(T) 



(exp(6u^)) 



-ll/s 11c 



-TU 



-s/11 



(102) 



. = -11/3, --^^exp(^-ATe)j= 0.718. 

In the region of intermediate temperatures, Eq. (102) agrees much better with the experimental data for the magne- 
tization of Fe and Ni [73] than the Bragg- Williams approximation. In distinction from the Bur ley extrapolation [74], 
our formula is simple and has a transparent physical background, taking into account both short-and-middle-range 
correlations, which contribute through the exponential function at low and intermediate temperatures, as well as 
long-range correlations dominating at the critical point. 



VIII. MAGNETIC PROPERTIES 



The ground state properties of the two-dimensional Heisenberg antiferromagnet can be considered by means of an 
expansion around its Ising limit in powers of the anisotropy parameter x, equal to zero for the Ising model and equal 
to one for the Heisenberg-limit [75,76]. Such expansions can bo generated separately for the most interesting cases 
of spin-1/2 and spin-1, allowing thus to treat them independently, not relying on the expansion in inverse powers of 
spin around a very distant case of classical spins. 

The expression for susceptibility x of the two-dimensional Heisenberg antiferromagnet was obtained in the following 
form [75,76]: 



2% ■ 



i + aix + a2X^ + asx^ + a^x"^ 



[x ^ 0), 



(103) 



ai 



a2 = 0.177083, as = -0.1898148, 04 = 0.191761, 05 = -0.196579, 



ae = 0.197934, 07 = -0.201447 



S = 



ai = -0.142857, 02 = 0.144643, 03 = -0.149916, 04 = 0.150672, 05 = -0.153095, 



ae = 0.153437, 



ar = -0.154932 



(5 = 1). 



According to general prescriptions of Section IV, we should, first, analyze the values of local multipliers mh{s), as 
s — > 00. 

Consider the case of 5 = 1. Since the values of the coefficients in the expansion (103) are slowly growing and 
oscillate with increasing number, the local multipliers will oscillate too. We conclude that in some cases we can 
continue the trajectory moving along the stable regions, while in other cases we have to move along unstable regions. 
In this situation we must rely on the aposteriori analysis. In order to choose the starting term in (103), let us compare 
the values of mfc(s), as s ^ 00. Then mi (00) = 1 + 02/01 = —0.143 and 1712(00) = 1 + 03/02 — —0.013. Since the latter 
number is smaller, we should start the renormalization procedure from the linear term, keeping the constant term 
untouched. The following sequence of exponential approximants can be readily written down, with the corresponding 
multipliers shown in brackets, calculated at r = 1 and x = 1: 



1 */ X 1 

2 X2{t) = g +aia;exp 



02 
ai' 



XT] = 0.146 {ml = 1), 



(104) 
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^ xHt) = ^ +aia;exp 



a2 


fa, W 


— xexp 


XT 


ai 


Ka2 J_ 



0.05 {ml = -0.025), 



(105) 



1 */ N 1 

2 Xiir) = - +aixexp 



02 






— a;exp 


— xexp 






Va2 





0.107 (mg = 0.944), 



(106) 



1 */ ^ 1 

2 XsW = g +aia;exp 



02 /as / a4 / as 
— a;exp I — a;exp — a;cxp — xt 
ai \a2 \a3 \a4 



= 0.075 



{ml = 0.225), 



(107) 



a2 (as ai a^ / ae 

— xexp I — xexp I — xexp I — xexp I — xt 
a\ \a2 Vfl3 \a5 



= 0.094 



{ml = 0.754), 



(108) 



2 = 8 + 



+aix exp 



a2 (as / a4 ( ae ( a-j 

— a; exp — a; exp — a; exp — a; exp — a; exp — xt 
ai Va2 Ka, ya^ Vas Vae 



0.083 



{ml = 0.394) . 



(109) 



We observe two subsequences, with odd and even numbers, probably embracing the correct result from below and 
above, respectively. We can suspect that they both define a stable quasifixed point corresponding to a focus. In order 
to locate it with maximal possible precision, we can impose a minimal difference condition on the points, belonging 
one to the "odd" and another to "even" subsequences, with the smallest absolute values of multipliers. On the other 
hand, it is interesting to compare the values that can be obtained from the minimal- difference condition imposed on 
two starting terms of the renormalized sequence with those obtained from two last terms. E.g., from the condition 
min^ |X3(''') ~ X2(''')l' obtain the value of the control parameter r = 0.56. Correspondingly, x^{t, 5 = 1)= 0.088. 
We obtain thus, 0.087 < %*(''": 5* = 1) < 0.088, with the lower bound following from the analysis of two last terms 
and of two approximants with minimal multipliers. Our estimate should be compared to the result of Ref. [76], 
0.095 ± 0.002, obtained by Pade-summation. 

The case of 5 = ^ may be treated by analogy with the case of 5 = 1. Following literally the same steps as above, 
we obtain 0.062 < x*(r, S = \) < 0.064, being close to the result of Pade-summation, 0.065 ± 0.003 [75]. 

The ground state energy E of the two-dimensional Heisenberg antiferromagnet may be presented in the form of an 
expansion [75,76]: 



2^ ~ -4 + a2a;^ + a4a;^ + aea;^ + a^x^ 



{x ^ 0), 



(110) 



a2 



04 = 0.0064, 

o 



ae = -2 X 0.00632628, ag = -2 x 0.0030085 ( 5 = - 



02 = -0.571428, 04 = -0.0504597, oe = -0.0144762, og = -0.00656238 {S = l). 



In the case of spin- 1, the values of the local multipliers mk{s) are increasing with increasing index k. In order to choose 
the starting term in (110), we compared the values of mi(oc) = 1 — a2/4 = 1.143 with 7712(00) = 1 + 04/02 = 1.088. 
I.e., an approximation cascade will have a more stable beginning if it starts from the second term in (110). Thus, we 
come to the exponential approximants 



2E2 = -4 + a2a;^ exp —x^t 



a2 



(111) 
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2ES = -4 + a2x'^ exp 



04 2 0,6 2 

— X exp — X T 

02 \a4 



(112) 



From the minimal-difference condition min,- |-E|(r) — E2{t)\ we find r = 1.568 and E^{t, S = 1) = —2.328, close to 
the result of Pade-summation, E = -2.327 ± 0.001 [76]. 

In the case of spin 1/2, the same approach is applicable, and we obtain r = 0.457 and AE^{t,S = 5) = —2.665. 
Pade-summation gives in this case, = —2.6785 ± 0.001 [75]. 

Pass now to considering magnetization M. In the case of spin-^ it is more convenient to present M in the form of 
a series in the parameter 5, related to x by the equation 



Then [75], 



l-5={l-x^f'\ 

2M ~ 1 + ai5 + a2(5^ + 03^^ + 04(5"* + asS^ {S 0), 



(113) 

(114) 



01 = --, 02=0.08, 0,3 = -0.009319, 04 = -0.4642, 05 = 0.08257. 

The local multipliers rrik behave quite irregularly, reflecting the behavior of the coefficients. In this situation we resort 
to an aposteriori analysis of the sequences of exponential approximants and corresponding aposteriori multipliers, 
whose values will be given in brackets. For the sequence including the constant term we observe a recurrent behavior, 
signaling an emergence of a limiting cycle; 



2Mi* = exp {ai6) = 0.641 {ml = 1), 



2M2* = exp 



oi^exp ( —5 

Oi 



= 0.69 (to2 = 0.737), 



(115) 
(116) 



2M3 = exp 



a2 



oi^exp — (5 exp — S 
oi \a2 



0-3 



= 0.685 {ml = 0.781), 



(117) 



2M4 = exp 



oioexp — oexp — dexp — 
oi V 02 Ka^ 



0.641 



{ml = 1), 



2M5* = exp 



Old exp — dexp — dexp — dexp( — 0) 
oi V 02 V03 04 



0.679 {ml = 0.808). 



(118) 
(119) 



We construct also a different sequence of exponential approximants, not including into the renormalization procedure 
the constant term: 



2M2* = l + oi(5exp ( ^(5 



0.629 



{ml = 1), 



(120) 



2M3 = 1 + oi^exp 



^2 . / 




— exp 




Oi 


V«2 j. 



2M4 = 1 + ai(5exp 



02 . / 03 / 04 
— exp — d exp — d 
oi \a2 \az 



2M^ = 1 + ai6 exp ( —5 exp ( —5 exp ( —S exp{—6) 
\ai \a2 \a3 ai 



0.621 {ml = 0.802), 
= 0.556 {ml = 5.149 x 10"*^), 
= 0.613 {ml = 0.722). 



(121) 



(122) 



(123) 



The last value agrees well with the result of Ref.[75], 2M = 0.605 ± 0.015. 
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In the case of spin-1, the following expansion in powers of 5 can be obtained [76]: 

M ~ 1 + 02(5^ + asS^ + 04(5* + 05^^ ^ 0), 



(124) 



02 = -0.326528, 03 = 0.326528, 04 = -0.73216, as = 1.3001056. 

The values of local multipliers m,k, in this case, suggest that inclusion into consideration of the last term destabilizes 
the trajectory. In this situation we again resort to the aposteriori analysis. The following sequence of exponential 
approximants can be readily written down, with the corresponding multipliers, calculated at r = 1, shown in brackets: 



M2* = exp (a25V) = 0.721 (mj = 1), 



(125) 



MS 



exp 



a2^^ exp ( —5t 
02 



0.887 (m; = 0.206) 



(126) 



Ml 



exp 



02^^ exp ( —5 exp ( —St 
a2 V«3 



0.746 {m*3 = 0.931), 



(127) 



Mg = exp 



a2d exp — oexp — oexp — or 
0-2 \az V"4 



0.848 {ml = 0.33). 



(128) 



We observe two subsequences, with odd and even numbers, probably embracing the correct result from below and 
above, respectively. We can suspect that they both define a stable quasifixed point, corresponding to a focus. In order 
to locate it with the maximal possible precision, we impose the minimal difference condition on two points, belonging 
to two different subsequences with the smallest values of multipliers, i.e.. min^ |M|(r) — M^{t)\. From this condition 
we determine r = 0.404 and M4 (r, 5= 1) = 0.804, in agreement with the estimate M = 0.81 ± 0.01 [76]. 

The expansions for the frequency moments, pi and p2, of the intensity of light scattering on the spin-pair excitations, 
are available for 5 = 1 [76], 



2pi 



ao + a2a; 



■ a^x 



(x^O), 



(129) 



ao = 14, a2 = 0.530612, 03 = -0.141138, 04 = 0.033837, as = 0.0171678; 

4p2 — ao + a2X^ + a^x^ + a^x'^ + a^x^ (a; ^ 0), 



(130) 



ao = 196, a2 = 17.469, 03 = -5.33453, 04 = 3.32992, as = 0.626767. 

Consider the case of p\. Comparing the local multiphers, as s ^ 00, mi (00) = 1 + 02/00 = 1.038 and 7712(00) = 
1 + 03/02 = 0.734, we conclude that the constant term should not be included into the renormalization procedure. 
Then, at r = 1 and x = 1, the following values of the exponential approximants can be calculated: 



2 {pI)2 = ao + ct2a;^ exp ( ^2;t 



7.203 {ml = 1), 



(131) 



/ 03 




— a; exp 




\a2 





7.215 {ml = 0.833), 



(132) 



2 (Pij4 = oo + 02a; exp ( —X exp ( —x exp ( —xt 



02 



as 



04 



7.222 {ml = 0.585). 



(133) 



We observe a sniootlily behaving sequence of multipliers. From the minimal difference condition minT- {{pD^ — ipl)^\, 
we find t = 1.59325 and (^1)4 = 7.221, in agreement with pi = 7.22 ± 0.02, quoted in Ref. [76]. Identical analysis 
leads to the value (^2)4 = 52.804 (r = 1.27), again in agreement with p2 = 53.0 ± 0.3 from Ref.[76]. The value of the 
parameter R = (-\/ P2 — Pi)/ Pi is equal to 0.113, close to 0.12 ± 0.3 from [76]. 
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IX. CRITICAL PHENOMENA 



A. Martinelli— Parisi e— Expansion 

Martinelli and Parisi suggested an interesting way to control the position-space renormalization group calculations 
[77], connecting the approximate Migdal-Kadanoff transformation with the exact theory by means of the control 
shift parameter e, equal to zero for Migdal Kadanoff approximation and equal to one for the exact renormalization 
transformation. This approach generates the expansions in powers of e, considered as a small parameter, around the 
Migdal-Kadanoff results. Finally, in order to reach the "exact" solution, one should set e = 1. The results can be 
further improved by imposing the condition on zero derivative of physical quantities at e = 1 [77,78]. 

The following expansion for the critical index v of the Ising model is available [77] for the square lattice: 

~ M;o + M;ie + u;2e^ (e 0), (134) 
wq = 0.687, wi = 1.14, W2 = -1.21. 

Directly for the index v, one can find 

~ ao + aie + 026^ (e ^ 0), (135) 

ao = 1.456, ai = -2.415, aa = 6.572. 

When only two starting terms from (135) are considered, the result is definitely wrong, v = —0.96, and with three 
terms we get v = 5.61. This shows how expansion (135) is bad. The renormalization procedure, not including into 
consideration the constant term, gives the exponential approximant 

j/*(e) = ao + aieexp ^^e^ = 1.297. (136) 

This result, is. probably, too large and it does change much, to the value 0.568, when the condition on zero derivative 
is imposed. Then, in order to extend the validity of the expansion (135), let us add to it a negative trial term — |o3| e^. 
The following renormalized expression can be written: 

2/3(6,03) = ao + oieexp I —e \ — -). (137) 

From the condition equivalent to (45), 

diy^{e,a3) 



de 



0, 



at e = 1, we find as = 4.268 and ^^{e = 1) = 0.992, in excellent agreement with the exact result, 1/ = 1. Our estimate 
is much better than the result of Pade-summation 1/ = 0.945, quoted in [77]. 
For the critical temperature Tc the following expansion was obtained [77] : 

~ 60 + 6ie + 626^ (e^O), (138) 
60 = 0.4359, bi = 0.024, 62 = -0.109. 

This gives for Tc 

Tc ~ ao + aie + a^e^ (e 0), (139) 

ao = 2.294, ai = -0.126, 02 = 0.581. 
Two exponential approximants can be written, 
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(^c)i = aoexp ( — er 



(140) 





ai 




ao exp 


— eexp 






ao 





(141) 



From the minimal difference condition, at e = 1, we find r = 0.278 and (T*)^ — 2.259. The result does not change 
much (to the value 2.248, r = 0.217) when the condition on zero derivative of (T*)^ is imposed. Let us add to the 
expansion (139) one more term — |a3| e^, so that the following approximant can be written: 



(^c)3 =«oexp 



ai a2 1 

— eexp — e , — 

ao \ai i + lMe 



an determine a-^ from the condition on zero derivative at e = 1. Then, \a^\ = 0.664 and (T* 



c J3 



(142) 



2.279. Fade 



approximants in this case give a close result, Tc = 2.275 [77]. This is to be compared with the exact Tc = 2.269. 
On a triangular lattice, the following expression for v was obtained [78] : 



2V- ^ a 

ao = 1.6786, 
In this case, the exponential approximant 



(2V^)*=, 



f aie + 


aae' (e ^ 0), 


ai 


0.5344, 02 = 




ai / 02 V 


a exp 


— eexp — e 




oo \ai J 



(143) 



-0.3952. 



(144) 



leads to v* = 1.035, which is a much better value than 1.161 obtained in Ref.[78] directly from (143). It is possible to 

improve our estimate performing the last step of the self similar bootstrap along the most stable available trajectory, 
with the stabilizer s corresponding to zero value of the local multiplier m = 1 + 02(1 + s)/ais. This yields 



Oo exp 



fli 
— ( 
ao 



. 01 J 



s = —- 



a2e 



ai + a2e 



= 2.e 



(145) 



Then, u* = 1.015 for r = 1. The derivative of z/* is equal to 0.035. At the point r = 1.152, the derivative goes to 
zero, and our estimate changes slightly to v* = 1.036. 

For the inverse critical temperature on a triangular lattice [78] , we have 



(146) 



60 = 0.3047, h = -0.0976, 62 = 0.0501. 
For the critical temperature we obtain 

Tc ~ ao + aie + a2e^ (e ^ 0), 



(147) 



ao = 3.282, ai = 1.051, 02 = -0.203. 

The following exponential approximant is favored from the viewpoint of local multipliers: 

T* = ao + aieexp (^^^t 

The derivative of (148) is quite large at r = 1, so we resort to the condition 

dT* 



(148) 



de 



0, 
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and at e = 1 we find r = —02/01 = 5.177, and 

T* =oo + aie"^ = 3.669, 

deviating from the exact value 3.642 with the percentage error 0.741%, while the result 3.888, following from expansion 

(146) and quoted in [78], gives the error 6.755%. 

For the first coefficient of the beta function /3i of the two-dimensional non-linear sigma-model the following 
expansion was obtained [79]: 

/3 

/3i~-^(ao+aie + a2e2) (e ^ 0), (149) 



1 3 69 



The exponential approximant 



V3 



0,2 



P*i = l^ao + aieexp (^-^e ) ) = -0.081 (150) 



is in good agreement with the exact result l^i = —j^ = —0.08 [79]. The derivative of (31 is small and equals to 0.006, 
so that we can safely stop at this point. 



B. Localization Length 

The critical exponent u describing the divergence of the localization length in the vicinity of the Anderson transition 
from the insulating to the conducting phase, 

/ ~ {E„ - £;)-^ 

as the energy E of an electron approaches the mobility edge E^, can be presented in the form of (2 + e)-expansion, 
or {d — 2)-expansion, where d is the dimensionality of space [80], 

1/ ~ i + 6e2 + (g ^ q), (151) 



q 27 
' = -4^(')' ^=T6^(^)- 

The result given by the starting two terms is wrong, v = —1.705, while it is known that u > 2/3 in the three- 
dimensional case (e = 1) [81]. The third term slightly improves the situation, but the result remains small, = 0.122. 
Renormalizing the starting two terms yields the exponential approximant 

z/*(e,T) = iexp(6eV) . (152) 

Let us impose the condition 

de " ' 

discussed above in Section IV. Then, the control function r(e) is found as follows: 




(153) 



At e = 1, wc obtain t = —0.123246, corresponding to the necessity to move backwards, from the wrong point of the 
approximation cascade. Our estimate for the critical index ly* = exp (^) = 1.39561 better agrees with the numerical 
result 1.3 [82] than the Pade-Borel estimate 0.730 quoted in Ref.[80]. 
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C. Amplitude Ratios 



Consider the 3d Ising model. Different amplitude ratios are available in the form of the Wilson e-expansion 
(e = 4 — (i) around the dimensionality four [83]. For the ratio /C~, related to the magnetic susceptibility in zero 
field, the following expression is available: 



C+ 



— ~ 2'>' [1 + aie + a2e^ + ase^] (e ^ 0), 



(154) 



1 



a2 



25 
108' 



1159 



11664 ' 

In order to improve the stability of the procedure, let us invert C+/C~ and study 

^) ~ 2-^ [1 + 6ie + 62e' + 63e'] (e ^ 0), 



A = 1.171953. 



(155) 



^^ = -2' 



b2 = 



54' 



63 = -7.5519 10" 



Prom the viewpoint of stability conditions, the following two approximants are well justified: 

-i\ * 



C+_ 

c- 



= 2-T 



1 + bieexp ( —CT 



(156) 



91 
c- 



1 + Oieexp I — eexp(— erj 
Oi 02 



(157) 



and from the minimal difference condition 



with the typical value of 7 = 1.24, we obtain r = 0.297. Correspondingly, (C^/C")* = 4.673. This value agrees well 
with the theoretical estimates and experimental data [83]. 

The ratio /A~ is related to the two-point correlation function at zero momentum, 



A+/A- ~ 2«-2 [l + e- [aal 



(e-0). 



(158) 



a2 



43 
54 



:A-C(3). 



Two renormalized expressions can be constructed, 



2°-^ [exp(eT)] 



(159) 



A+ 



2"-2[exp(eexp(-|a2|eT))] 



(160) 



both exponential approximants being justified from the viewpoint of stability conditions. Prom the condition 
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with the typical value of a = 0.11, we find that r = 0.669 and = 0.527, well agreeing with the data of Table 

5 from Rcf. [83]. 

For the amplitude ratios Rc and (see [83]), the following expansions are available: 



^ i2-2/3-ie[i + aie-]a2le^] (e ^ 0), 



(161) 



ai 



17 
27' 



989 4^,,, 2^ 
"^=2916-9^(^^-3^' 



and 



R ~ 3('5-3)/227+(l-5)/2 



(e-0). 



(162) 



The self-similar exponential approximants are 



(i?c); = -2-2/3-ie[exp(aieT)], 



(163) 



exp I aieexp 



ai 



-er 



(164) 



From the minimal difference condition min^ |(-Rc)2 ~ (-^c)i|) with the typical value /? = 0.325, we obtain r = 0.633 
and {Rc)*2 = 0.053 agreeing well with the results quoted in Ref.[83]. 
The renormalized expression for R^^ writes 



R* = 3('5-3)/22T+(l-'5)/2 



(165) 



We take S = 4.825, corresponding, by virtue of the scaling relation, to a reasonable value for the critical index r] = 0.03. 
Then i?* = 1.675, close to various estimates presented in [83]. 

Finally, we consider the quantity |2:o|, the universal rescaled spontaneous magnetization, 



^0 



V32^ [l+aie + aze^- jaaje^] (e ^ 0), 



(166) 



1 73 1 \ 7 yro\ 5581 

"^ = 4' ""'=864' "^=24^ "36^(^^ + 93312- 

The self-similar exponential approximants, justified from the viewpoint of an aposteriori stability analysis, are 



\zq 



; = ^/32'5 



exp [ aieexp [ — er 



(167) 



1^0 Is = V32^ 



02 

exp I aieexp I — eexp 



(168) 



and from the minimal difference condition min^ ll-^ols — |^ol2| find r = 0.633 and |2;o|3 = 2.913 agreeing well with 
the result of Pade-summation 2.87 ± 0.06 [83]. 



D. Critical Indices from the Wilson e— Expansion 



Critical indices are usually obtained from the Wilson e— expansion [84,85] using some kind of a resummation 

procedure. As a rule, for different values of n, standing for the number of the order parameter components, one 
obtains some values not related to each other analytically. We obtain below analytical renormalized expressions for 
the critical indices u and r/, valid for arbitrary n. 
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For the critical index 77, the expansion is available up to the fifth order term in e [85]. For convenience, we reproduce 
the higher-order coefficients in Appendix, 



r] ~ a2(n)e^ + a3(n)e^ + ai{n)e'^ + a5(n)e^ (e 0), 



(169) 



02 (n) = 



n + 2 



n + 2 



■(272 + 56n-n2), 



2(n + 8)2' 8(n + 8)4' 

The following exponential approximants are justified from the viewpoint of the stability conditions for local multipliers: 



r]l{e,n,T) = a2(n)e^ + a3(n)e^ exp 



04(71) 
03 (n) 



(170) 



r?5(e,n,T) = 02(71) + 03(11) exp 



04(71) / 05(71) 

eexp I — r^er 



03(71) \a4,{n) 

From the minimal difference condition min^ IVai^j 't) ~ Vli^^ it-jt)] , equivalent in this case to the equation 

'05(7^) 
r = exp I — r^er 



(171) 



04(n) 

wc can find the control function r = ''"(ti) at e = 1. From formula (171), wc can estimate the critical index 
7^5(1, 71, T(n)) = r]*{n). The results of calculations are presented in the Table. They agree well with the majority 
of different kinds of estimates available for the critical index rj. 

For the critical index u, up to the fifth order in e (see Appendix), one has [85] 



u-^ ~ bo{n) + bi{n}e + b2{n)e^ + b3{n)e^ + bi{n)e'^ + b5{n)€^ (e ^ 0), 



(172) 



6o(n) = 2, bi{n) = - 



n + 8' 



62 (n) = - 



(n + 2)(137i + 44) 
2(71 + 8)3 



The following exponential approximants, preserving the correct limits v = ^ a,t n = —2 [86- 
[89] , can be written, giving at r = 1 and e = 1 the following results 



(v ^)* = bo{n) + 6i(n)eexp 



&2(n) 
61 (n)' 



and = 1, as n — > 00 



(173) 



i/*(n = 0) = 0.607, z/2(n = 1) =0.655, 2/3 (n = 2) = 0.698, 1/2(71 = 3) = 0.736. 



bo(ji) + &i(n)eexp 



eexp I . , : er 



61 (n) 



b2{n) 



(174) 



i/*(n = 0) = 0.579, 1/3(71 = 1) = 0.616, 7/3 (71 = 2) = 0.651, 1/3(71 
(z/"^)* = 60(71) + 6i(7i)eexp 



&2(7l) A3 (71) A4(7l) 

eexp ———eexp | . , , er 



61(71) V^2(7l) 



63(71) 



3) = 0.683. 



(175) 



i/*(n = 0) = 0.603, z/|(7i = 1) = 0.649, z/|(7i = 2) = 0.692, 1/^(71 = 3) = 0.729. 
(z/~^)* = 60(71) + 6i(7i)eexp 



62(71) ^63(71) /64(7l) /65(7l) 

eexp I eexp ( e exp ( . , , er 



61(71) V^2(71 



63(71) 



64(71) 



(176) 



z/*(n = 0) =0.58, 1/5(71 = 1) = 0.618, f5(7i = 2) = 0.654, 2/^(71 = 3) = 0.688. 
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Two sequences, with odd and even numbers, are clearly seen. The last two approximants give the closest values. 



From the minimal difference condition mini- 



which simply reduces to the equation 



'hz(n) 



one can easily find the control function r = T(n), and, finally, calculate the critical index, v*{n) — (^{y 

For the critical exponent w(e) = 2B^, calculated in the renormalized infrared-stable fixed point, the following 
expansion is available [85]: 



uj{e) e + C2{n)e^ + cs{n)e^ (e ^ 0), 



(177) 



9n + 42 
(n + 8) 



4(n-|-8)'* 



[33n3 + 53871^ + 4288n + 9658 + C(3) (n + 8) 96(5n + 22)] . 



Wc limit here the discussion by the third-order terms, since the higher order terms grow impetuously, jeopardizing 
the fulfillment of the stability conditions. The following exponential approximants are available being well justified 
from the viewpoint of stability, 



w2(e, n,T) = eexp (c2(n)eT) 



(178) 



uj'^{e, n, t) = eexp 
and from the minimal difference condition 



C2(n)eexp — ^er 
C2 n) 



(179) 



T = exp — r^er 
\C2{n} 

one can find the control function r = r(n) at e = 1. We obtain the following values of the index in the physically 
interesting region: 

a;*(n = 0) = 0.788, a;*(n = 1) = 0.788, a;*(n = 2) = 0.791, a;*(n = 3) = 0.794. 

The limiting values, ui*{n = —2) « 0.8, u}*{n — > oo) ~ 1, sound reasonable. The dynamical critical index 6 = ivu can 

be estimated using (176) and (179). The results for the critical indices i] and i/, presented in the Table, well agree 
with other theoretical estimates and available experimental data [90-92]. Let us stress again, that at n = —2 and 
n — > 00, we obtain the exact results. 



E. Critical Indices from the Field Theory Expansion 

Field-theory approach in the theory of critical phenomena is, usually, very accurate [92,93]. In Ref [13] we analyzed 
the expansions in powers of the interaction constant g (g-expansion) for some critical indices from the viewpoint of 
the limiting cases n = —2, n ^ oo, and found by a direct inspection of the expressions for ij and 7 from [92], that 
the n — > 00 limit, corresponding to the spherical model [89], is obeyed rigorously if 51 = 1, i.e. r] = 0, 7 = 2, and the 
n = —2 limit, corresponding to the gaussian polymer [86-88], is obeyed with a very high accuracy for arbitrary g, i.e. 
ry « 0, 7 « 1. 

The standard approach [93] uses, for computing the renormalized infrared stable fixed point g* of the beta-function 

Bi{g), a complicated Borel summation technique. Then critical indices are calculated as 7(3*) and ri{g*). The results, 
thus, depend on the way in which the position of the fixed point is determined, although, different approaches give 
the results very close to each other [92]. We suggest below a simple way to minimize an uncertainty related to the 



27 



position of the fixed point. Let us use below two different approaches to the determination of g*. The first approach 
was suggested in Ref.[f3]. It is based on three starting terms from the g-expansion from Ref.[92] (see Appendix): 

~ -5 + £,2 + as (n)ff 3 + 04(71)5'^ + 05(71)3^ + 06(71)/ (180) 
, , 6.07407408 71 + 28.14814815 

«3W = (^TTsp ' ••• 

For the velocity field we have 



«3(/) = -f (i(i + 7rT47y 



which is to be substituted into the evolution integral 

df 



T. (181) 

9+9^ «3(/) 

The root gl = gl{n,T) of the equation Bl{g,n,T = 1) = is obtained numerically, as a function = gtin)- The 
following values were obtained in the physically important cases: 

.gi(r2 = 0) = 1.59, .gi(n = 1) = 1.559, .gi(n = 2) = 1.524, gl{n = 3) = lA91. 

Such a decreasing, with n, dependence is characteristic to the majority of related studies. We observed also, that at 
n = —2, gl = 1.599 and at n — > 00, gl = 1. The dependence of g*{n) in the interval n G (—2,0) is nonmonotonous, 
a maximum is reached at 7i = —1, where gl = 1.61. We also constructed a different beta-function -82(5), 

^2(5. t) = -5 exp (-03(71)9^) + 5^ exp (04(71)31^) + 05(71)5^ exp ( ^^^4~t5M > (182) 



^05(71) 

leading to the following increasing with n values of the fixed point g2{n, r) at r = 1, 

5^(71 = 0) = 1.163, 52*(ri = l) = 1-208, 5^(71 = 2) = 1.249, 5^(n = 3) = 1.285, 

possessing a maximum at n = 7. Such an increasing dependence better corresponds to the known decrease of the 
critical temperature Tc , from the Ising (n = 1) to Heisenberg (n = 3) model [72], since Tc ^ 9^^ [94]. We attempt 
now to minimize uncertainty, connected to the way of determining the fixed point. Imposing the minimal difference 
condition 

min|5f2(7i,r) - gl{n,T) \ , 

T 

we obtain the control function r = T(n), and the following values for the optimized zero of the beta-function 

£,*(n = 0) = 1.311, 5*(n = 1) = 1.334, g*(7i = 2) = 1.352, 5*(n = 3) = 1.365. 

These values increase till n = 5, when g*{n = 5) = 1.372, and then decrease till 5* = 1, as 71 — > 00. 

For the critical index 77, we keep all the terms available, up to the sixth order in powers of g [92] (see Appendix): 

r]{g) ~ 62(^)5' + &3(n)5' + bi{n)g'' + 65(^)5' + h{n)g^ {9 ^ 0), (183) 
, , , 0.2962962963(71+2) , , , 0.024684001471^ + 0.24684001471 + 0.3949440224 

^^(") = ^7T8? ' ^^(") = ^^Tsp 

From the aposteriori stability analysis we select the fifth order approximant corresponding to the smallest multiplier: 



nlig^n) = 62(^)5^ exp 



"5exp 7---gexp 7— ^ff 



b2{n) \h{n) \b4{n) 



(184) 
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so that 



r?5*(ff*,n = 0) = 0.027, r?5*(ff*, n = 1) = 0.034, r]*^{g* ,n = 2) = 0.039, r]*^{g* ,n = 3) = OM. 



These values agree well with those quoted in Refs. [91,92]. 

For the critical index 7 we write down all the terms available, in powers of g [92] (see Appendix), 

7-1 ~ 1 + ci{n)g + C2(n)ff2 + 03(71)5^ + C4(n)ff^ + C5{n)g^ + ce{n)g^ (5^0), 



(185) 



ci(n) 



n + 2 
'2{n + 8)'' 



C2{n) 



n + 2 
(n + 8)2' 



The third order exponential approximant is selected by an aposteriori analysis, since it corresponds to the smallest 
value of the multiplier, giving the following expression for the critical index: 



(7 ^)* = l + ci(n)5exp 



C2(n) / C3{n) 
—r^g exp — 
ci(n) \C2{n) 



(186) 



so that 

7*(5*,n = 0) = 1.164, 7*(5*,n = l) = 1.243, 7*(ff*,n = 2) = 1.319, 7*(ff*, n = 3) = 1.39, 
which again agrees well with the data of Refs. [90-92]. 



X. DISCUSSION 



In this paper we have developed an analytical approach for summing divergent series with arbitrary noninteger 
as wcill as integer powers. This approach is based on the novel notion of the self similar exponential approximants. 
By a number of examples we show that the developed method is general and accurate. In addition, the exponential 
approximants have a simple analytical structure, even for quite large number of perturbative terms used, when one 
usually has to resort to numerical techniques. Because of its analytical nature, our method permits one to accomplish 
direct analysis of resulting formulas with respect to the variation of physical parameters. All examples we have 
analyzed are related to important physical phenomena, such as critical and crossover phenomena. We demonstrate 
that the method of self-similar exponential approximants provides an effective general tool for treating many different 
problems of statistical physics. 

In conclusion, it is worth touching the following question. Assume that we are given a truncated series (1). Then, 
would it be possible, looking at the given series, to get some heuristic arguments, when the self-similar approximants 
in the form of nested exponentials should yield good results. The answer to this question is: Yes, we can made a 
preliminary estimates whether the nested exponentials would work well. Such a first-glance investigation can be done 
by analysing the apriori multiplier (42). Since the case of the nested exponentials corresponds to s ^ 00, then from 
(48), we have 

k 

lim mk(x, s) = y — a;""""". (187) 

From here, we immediately notice what would be the favourable cases for the better stability of the procedure, which 
is related to the condition of the minimal multiplier modulus, |rnfe(a;, oo)|. Such favorable cases include: (i) When 
a„ decreases as n increases, so that |a„/ao| < 1 and a„/ao — > 0, as n ^ 00. If |a:;"'*~"''| increases with n, then, to 
compensate this increase, |a„/oo| must decrease sufficiently fast, (ii) When with increasing n, \x°'"~°'° \ decreases. This 
kind of situation occurs, e.g., in the strong-coupling limit of many quantum problems, when a„ < 0, and x""""" — > 0, 
as n ^ 00. If the decrease of is sufiiciently fast, then |a„/ao| may even grow with n. (iii) When (1) is an 

alternating series, that is, the coefficients a„ change their signs with changing n. In this case, even if neither ja^i/aol 
nor |a;""~"° | decrease, but nevertheless, because of alternating signs of a„, the value \mk{x, oo)| may be small. 

The worst case, as is seen from (187), would be when all coefficients a„ are of the same sign, a„ increases with n, 
so that Un/ao > 1, and in addition, when x""~°'° > 1. To explain why this case is really the worst, consider a simple 
illustration for series (1) truncated at second order, 

P2{x) = ao + aix + a2x'^. 
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The related apriori multiplier (187) is 



m2{x,oo) = l + —x+—x^. (188) 

(12 do 

Starting with P2{x), let us construct the simplest nested exponentials, not invoking control functions, 

f*{x) = ao exp ( —X ] , /2 (a;) = ao exp ( —xexp ( —x ) ] . (189) 

The sequence {/^(a;)} has to satisiy the stability condition (50), with the aposteriori multipliers (49). In our case, 
m*(a;) = 1 and 

™«.)=(l.%)exp(H,)|M. ,1,0, 

Assume now that a„+i > a„ > and a; > 0. Then /*_|_;^(a:) > fn{x)- In such a case, multiplier (190) is more than 
unity. Hence, the procedure is locally mistable, and we cannot trust to approximants (189). The first of them, i.e. 
fi{x), can yet give a reasonable estimate, but the second. /I (a;), is certainly luitrustable. 

In this way, a quick glance at the apriori multiplier (187) gives us a feeling of whether the self-similar exponentials 
would produce good results. But the final conclusion of whether we did have managed to construct a convergent 
sequence of nested exponential approximants is to be based on the stability analysis of the renormalized multipliers 
(49). The relation between the latter and the apriori multipliers (187) is not direct, as can be seen even from the 
simplest case resulting in (188) and (190). In general, it would be more correct to say that there is no direct relation 
between these two types of multipliers. Therefore, there is no necessity of requiring that the apriori multipliers (187) 
be compulsory less than unity in their absolute values. It is sufficient to require that the final multipliers (49) satisfy 
the stability condition (50). 

Another question that may arise is as follows. Suppose, we have met with the worst case, when we are not able to 
construct a convergent sequence of self-similar exponential approximants. How then should we proceed in order to 
define an effective limit of a divergent sequence {pk{x)}7 In such a case, there are several possibilities. First of all, 
since it is always ml{x) = 1, as is clear from definition (49), thence fi{x) may serve as an estimate for an effective 
limit of the sequence {pk{x)}. 

Second, recall that the nested exponentials are only one of admissible forms of self-similar approximants, corre- 
sponding to a particular case, when the power in the algebraic transformation (2) is assumed to tend to infinity. If 
the latter assumption is waved aside, we return to the radical form (10). Then, we have to define control functions 
s = Sk{x) and r = Tk{x) so that the sequence {p^(a;, s^, Tfe)} be convergent. Equations defining these control func- 
tions, as always in the framework of the self-similar approximation theory [7-14], follow from the requirement for the 
corresponding approximation cascade to have a stable fixed point; the stability of a fixed point being in one-to-one 
correspondence with the existence of a limit which the sequence of approximations converges to [10,11,95]. Developing 
this procedure, we come to a convergent sequence {/^(a;)} of self-similar approximants f^{x) = pl{x, Sk{x),Tk{x)) 
which can give very accurate approximations for the sought function, but will not have such a nice structure as that 
of nested exponentials. 

The self similar exponential approximants, in addition to having a nice and convenient mathematical structure, 
evidently illustrate by their form the idea of self-similarity of approximations [7-9]. Thus, if we introduce a function 

G{x, y) = xe^ 

and use the notation 

Xi = biX^*; xo = aox°"', 
then the self-similar approximant (18) can be written as 

Fk{x) = G{xo,G{xi,. . . G{xk-i,Xk)) ■ . .). 

Finally, one more possibility of treating divergent series, when the direct construction of nested exponentials does 
not work, could be either by resorting to a change of variables or by invoking a transformation of the given series, 
so that the following application of the self-similar approximation theory would result in a convergent sequence of 
nested exponentials. What kind of a transformation or a change of variables is appropriate can again be decided by 
means of the apriori multipliers (187). 
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Appendix 

The higher-order coefficients in the e-expansion for the critical index rj are [85] 

, , 1 n+2 

X [5n^ + 230n^ - 1124n^ - 17920n - 46144 + C (3) (n + 8) 384 (5n + 22)] , 

1 n + 2 

X [{l3nP + 946n^ + 27620n^ + 121472n^ - 262528^2 _ 2912768n - 5655552) ■ 
-C (3) (n + 8) 16(n^ + lOn^ + 1220n^ - 1136n^ - 68672n - 171264)+ 

+C(4) (n + 8)^ 1152 (5n + 22) - C(5) (n + 8)^ 5120(271^ + 55n + 186) 
For the critical index v one has 

^3(") = „|""'"!L [Sn^ - 452n2 - 2672n - 5312 + C(3)(n + 8)96(5n + 22)1 , 
8{n + 8) 

64 (n) = ^2{ri^_^8y ^^"^ ^^^""^ ~ 1290071^ - 8155271^ - 219968n - 357120+ 

+C(3)(n + 8)16(3n* - 194n^ + 148n^ + 9472n + 19488)+ 
+C(4)(n + 8)^288(571 + 22) - C(5)(n + 8)21280(271^ + 557i + 186)], 

^^^"^= 128(n + 8)^ " 

X [371'' - 119871^' - 2748471^ - 105534471^^ - 524211271^ - 525670471^ + 699904071 - 626688 - 

-C(3)(7i + 8)16(137z'^ - 31071^ + 1900471'' + 10240071^ - 38153671^ - 279257671 - 4240640)- 

-C^(3)(n + 8)^1024(271* + 18n^ + 98171^ + 6994n + 11688)+ 

+C(4)(7i + 8)^48(371"^ - 19471^ + 14871^ + 947271 + 19488)+ 

+C(5)(7i + 8)^256(15571^' + 302671^ + 98971^ - 6601871 - 130608)- 

-C(6)(n + 8)^6400(271^ + 557i + 186) + C(7)(n + 8)^56448(1471^ + 18971 + 526)] . 
The higher-order coefficients of the field-theory expansion [92] for the /3-function are 

04(71) = - — ^—^(1.3489427671^ + 54.9403769871+ 199.6404170), 

(71 + 8)-^ 
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a^in) = — 5— (-0.15564589n^ + 35.82020378n^ + 602.5212305n+ 1832.206732), 

X (0.0512361871"* + 3.23787620n3 + 668.554336871^ + 7819.564764n+ 20770.17697) . 
For the critical index r], 

b4{n) = - — ^(-0.0042985626n^ + 0.6679859202n^ + 4.609221057n + 6.512109933, 

65(71) = -- — 3— T (0.006550922271''- 0.1324510614n^+ 1.89113928271^+ 15.1880934071+ 
(n + 8)5 ^ 

+21.64720643), 

^^{n) = 7 — 5—^ (-0.00554892071^ - 0.020399448577^' + 3.05403098777^ + 64.0774465677^+ 
[n + 8)° ^ 

+300.720893377 + 369.7130739) . 

For the critical index 7, 

C3(n) = -- — ^(0.8795588926n^ + 6.485476868n+ 9.452718166), 

C4(n) = - — 5— ^(-0.128332104377^ + 7.96674070377^ + 51.8442129877+ 70.79480631), 

05(77) = -- — 3—^ (0.049096605877^ + 4.28815249371^ + 108.36182197i2 + 537.813610571+ 
(71 + 8)^ ^ 

+675.6996077) , 

C6(n) = - — 3_ (-0.025926794571^ - 1.61862784371^' + 85.5456974671^ + 1538.8182357iV 

+6653.956526n+ 7862.074086) . 
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Table Caption 

Critical indices, for the models with different numbers of components n, calculated by using self-similar exponential 
approximants obtained from the e-expansion, and compared to the results listed in literature. The indices u, t], and 
6 are calculated directly, and other indices are found from the known scaling relations. 
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